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© Method of shading a graphics Image. 



© A method* and apparatus for tiling a display area 
defined by lines joining vertices with a surface nor- 
mal vector associated with each vertex whereby the 
display area is subdivided into sub-display areas by 
calculating additional vertices and surface norma! 
vectors by interpolation and rendering a given sub- 
display area by calculating intensity values at its 
vertices and tiling its area by linear interpolation of 
the calculated vertex intensity values. 
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METHOD OF SHADING A GRAPHICS IMAGE 



Background of the Invention 

This invention relates to a method of shading a 
graphics image. 

In a computer graphics system, an image is 
displayed on the screen of a cathode-ray-tube 
(CRT) or other display device. Generally, the elec- 
tron beam of the CRT Is deflected under control of 
a digital signal having an X (horizontal) component 
and a Y (vertical) component, so that pixels on a 
rectangular array of grid points are addressed row- 
by-row in accordance with a raster pattern. Data 
representing the image is generated in a picture 
processor and is down-loaded from the picture 
processor to a frame buffer. The frame buffer has 
memory locations which correspond on a one-to- 
one basis with the grid points of the CRT display. 
In the following description, it will be assumed for 
the sake of clarity that the memory locations of the 
frame buffer are physically arranged in a rectangu- 
lar array which corresponds to the array of grid 
points of the display. 

The digital signal used to control deflection of 
the CRTs electron beam is also used to control 
addressing of the frame buffer, and the color with 
which the pixel at a grid point is displayed depends 
on the contents of the corresponding memory loca- 
tion. The frame buffer thus contains a virtual image 
of the display that is provided by the CRT. The 
data provided to the frame buffer by the picture 
processor may represent an image composed of 
colored areas. The process of loading data values 
into the frame buffer memory to represent a col- 
ored area is called tiling. In the following descrip- 
tion, references are made in connection with a 
tiling operation to the tiling of specified grid points 
and the tiling of triangles. However, this is intended 
merely as an abbreviated way of referring to the 
tiling of the memory locations corresponding to the 
specified grid points and to the tiling of the mem- 
ory locations corresponding to the grid points that 
lie within the triangles. 

A computer graphics system may be used to 
depict a three-dimensional curved surface by re- 
solving the surface into triangular facets and dis- 
playing the projection of those facets into a two- 
dimensional display plane. Each triangular facet 
forms a triangle or a line segment when projected 
into the display plane. The case in which the pro- 
jection is a line segment will not be considered 
further herein. Information defining the triangular 
facets in three-dimensional world space is specified 
to the picture processor, which utilizes that in- 
formation to generate information defining the cor- 



responding display plane triangles. This information 
is made available by the picture processor to a 
rendering engine, which is used to load shading 
intensity values I into the frame buffer. Generally, 
s the triangular facets are treated sequentially by the 
picture processor. 

in order to provide an accurate depiction, it is 
necessary to take account of the fact that when a 
three-dimensional object is examined, surfaces of 
/o the object that are on the far side of the object 
from the viewer cannot be seen, and therefore data 
representing these surfaces should not be loaded 
into the frame buffer. In a graphics system, this is 
accomplished by hidden surface removal. Hidden 
is surface removal is discussed in Newman and 
Sproull, "Principles of Interactive Computer Graph- 
ics" (McGraw-Hill Book Company, 1979), Second 
Edition, 369. 

Data loaded into the frame buffer memory for 
20 tiling of a given grid point represents the color and 
intensity of light reflected towards the viewer by 
the corresponding point of the surface being de- 
picted. In the present discussion, this data is re- 
ferred to as an intensity value I, but may actually 
25 be several values representing intensities of sepa- 
rate color components, e.g. intensities of red, green 
and blue. The three-dimensional appearance of a 
surface will generally be depicted with an en- 
hanced degree of realism if the intensity value I for 
30 a grid point takes account of the position and 
orientation of the surface element that is depicted 
by that grid point. This may be accomplished by 
calculating a surface normal vector at the surface 
element and applying a lighting model to derive an 
35 intensity value I for a grid point representing that 
surface element. It should be appreciated that use 
of a lighting model to calculate a grid point inten- 
sity value is computationally expensive as complex 
calculations are required to account for light sour- 
40 ces and surface reflectivity. 

The intensity value I that is loaded into the 
frame buffer for a given grid point depends on a 
shading model that is applied to the displayed 
surface. Three major shading models are described 
45 in J. D. Foley and A. Van Dam, "Fundamentals of 
Interactive Computer Graphics", 580-584, (1984). 
The three models are referred to as constant shad- 
ing, intensity Interpolation, or Gouraud, shading and 
normal-vector interpolation, or Phong, shading. 
so Constant shading uses a single intensity value for a 
given display area and is much less satisfactory for 
depicting a curved surface than either Gouraud 
shading or Phong shading. 

In Gouraud shading of a display plane triangle 
that represents a triangular facet, a vertex surface 
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normal is given for each vertex of the facet and a 
vertex intensity is found for each corresponding 
display plane vertex by applying a selected lighting 
model to the corresponding vertex normal. The 
calculated intensity values provide a basis for 
shading the display plane triangle by linear inter- 
polation of these intensities along each edge of the 
triangle and between edges along each scan line. 
Thus, to Gouraud shade a triangle the selected 
lighting model is accessed once for each vertex. In 
Phong shading of a triangle, the vertex surface 
normals are given. End point normals are cal- 
culated for the points of intersection of each scan 
line with the edges of the triangle, by interpolation 
between the vertex normals at opposite ends of the 
edges that are intersected by the scan line. For 
each pixel along a scan line, a surface normal 
vector is calculated by interpolation between the 
end point vectors for that scan line. A surface 
normal vector is associated with every pixel and a 
selected lighting model is applied to the surface 
normal vector for each and every pixel to calculate 
its shading intensity. Because the lighting model is 
computed for every grid point within the triangle, 
Phong shading generally requires much more time 
to execute than Gouraud shading. However, Phong 
shading provides a more realistic depiction of an 
object because the position and orientation of each 
surface element corresponding to each grid point is 
taken into account. 



Summary of the Invention 



A preferred embodiment of the present inven- 
tion is a method of loading a memory device 
having a plurality of addressable memory locations, 
each of which can be defined uniquely by an 
address having an X component and a Y compo- 
nent, which locations correspond respectively to 
grid points in a rectangular array, and wherein it is 
desired to load or tile memory locations lying with- 
in a display area defined by at least three vertices 
(XO f YO), (X1. Y1) and (X2,Y2), and their associated 
surface normal vectors. The method comprises de- 
fining points (X3.Y3) and (X4,Y4) on the lines that 
join the point (XO.YO) with the points (X1,Y1) and 
(X2.Y2) respectively, such that a sub-area is de- 
fined by the points (XO.YO), (X3.Y3) and (X4,Y4). 
Surface normal vectors for points (X3.Y3) and 
(X4.Y4) are derived by interpolation of the surface 
normal vectors associated with points (XO.YO). 
(X1,Y1) and (X2.Y2). intensity values IO, 13 and 14 
for the points (XO.YO), (X3.Y3) and (X4.Y4) are 
calculated by application of a lighting model to 
their associated surface normal vectors, and mem- 
ory locations corresponding to points within the 



sub-area are loaded with intensity values that are 
derived from IO, 13 and 14. 

In one aspect of the preferred embodiment, a 
"parent" triangle presented for tiling models a por- 

s tion of a curved surface, the parent triangle being 
defined by three vertices and three corresponding 
surface normal vectors each taken as being normal 
to the curved surface. Congruent sub- triangles are 
defined by bisecting each side of the parent tri- 

io angle and surface normal vectors are derived for 
each vertex of each sub-triangle by interpolation of 
the surface normal vectors of the parent triangle. 
Each sub-triangle is tested to determine suitability 
for Gouraud shading. Sub-triangles not considered 

is suitable for Gouraud shading are further subdivided 
and the new sub-triangles are tested to determine 
suitability for Gouraud shading. Sub-triangles con- 
sidered suitable for Gouraud shading are tiled by 
first applying the vertex surface normal vectors to a 

20 lighting model to obtain vertex intensity values and 
interpolating these intensity values across the inte- 
rior of the sub-triangle. 

25 Description of Drawings 

For a better understanding of the invention, and 
to show how the same may be carried into effect, 
30 . reference will now be made, by way of example, to 
the accompanying drawings in which, 

FIG. 1 is a block diagram of a computer 
graphics system, 

FIG. 2 illustrates a rectangular array of grid 
36 points and a triangle having its vertices at grid 
points, 

FIG. 3 illustrates the triangle of FIG. 2, di- 
vided into sub-triangles. 

40 

Detailed Description 

The computer graphics system illustrated in 
45 FIG. 1 comprises a picture processor 1 which re- 
ceives a display list and generates information re- 
presenting triangles and applies information defin- 
ing these triangles sequentially to an XY address 
generator 2 and an I generator 6. XY address 
so generator 2 is described and claimed in EP-A2- 
0314367 and in EP-A2-03 14368, the disclosures of 
which are hereby incorporated by reference herein. 
I generator 6 is described in EP-A2-0314368. XY 
address generator 2 executes a modified Sressen- 
55 ham algorithm to select grid points within a display 
plane triangle by providing an X address in register 
3 and a Y address in register 4 for addressing a 
frame buffer 5. I generator 8 performs Gouraud 
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for each pair of coordinates (X,Y) provided by XY 
address generator 2. In Gouraud shading, it is 
assumed that the value of i varies linearly in both 
the vertical direction and the horizontal direction. 
Therefore, the intensity values 10, 11 and 12 for the 
vertices VO, V1 and V2 satisfy a plane equation in 
X and Y from which a horizontal intensity increment 
and vertical intensity increment may be calculated. 
Initially, picture processor 1 loads the horizontal 
and vertical intensity increment values in I gener- 
ator 6 and calculates an initial intensity value I for 
determining the intensity of the first grid point to be 
tiled, and this value is also loaded into I generator 
6. 

After initialization, XY address generator 2 and 
I generator 6 enter the operating state. In this state 
generators 2 and 6 operate synchronously to tile 
the triangle according to the modified Bresenham 
algorithm and the Gouraud shading model. XY ad- 
dress generator 2 calculates the X coordinate for 
the first pixel to be tiled on each row and the 
number Hrz of pixels to be tiled on that row, and 
also determines the direction (left to right or right to 
left) in which tiling is to take place. I generator 6 
provides an I value for each grid point addressed 
by generator 2. It will be noted from the co-pending 
Patent Application Serial No- 07/113,031 that no 
value of Hrz is calculated at the bottom grid row of 
the triangle and therefore the bottom grid row is 
not tiled. However, the top row is tiled. Similarly, 
the grid point at the floor of the exact intersection 
of a grid row with the right siding (V0V1V2 in the 
case of FIG. 2) and the grid point at one plus the 
floor of the exact intersection of the grid row with 
the left siding and all the grid points lying between 
the above two grid points on the given grid row are 
tiled. Accordingly, when triangles abut. e.g. in a 
tessellation, a grid point at the boundary between 
two adjacent triangles is treated as being in one, 
and only one, triangle. This ensures that there are 
no gaps between triangles and no overlaps of tri- 
angles. 

Gouraud shading the triangle V0V1V2 in this 
manner would be satisfactory if the portion of the 
curved surface represented by the triangle were 
substantially planar. However, this is not generally 
a valid assumption, and therefore Gouraud shading 
in this manner would not normally be satisfactory in 
relation to that obtained by Phong shading. On the 
other hand, Phong shading is computationally ex- 
pensive and often not necessary for a reasonably 
accurate rendition of a curved surface. 

In order to avoid this limitation on the useful- 
ness of Gouraud shading and avoid the computa- 
tional costs associated with Phong shading, a com- 
bination of Phong and Gouraud shading techniques 
is used. The display plane triangle is divided into 
four congruent sub-triangies (FIG. 3) by bisecting 



each edge of the triangle and connecting the bisec- 
tion points to form four sub-triangles. For each sub- 
triangle three vertex data structures are established 
with each structure containing the same param- 

s eters as the vertex data structures for vertices of 
the parent triangle. The positional coordinate values 
and surface normal vector components for the new 
vertices are obtained by interpolation. For example, 
the bisection points of the edges V0V1, V1V2 and 

70 V0V2, which may be designated V01 f V12 and 
V02, respectively have coordinate values interpo- 
lated from the coordinate values for points V0 f V1 
and V2, and have surface normal vector compo- 
nents interpolated from the surface normal vector 

is components for points VO, V1 and V2. A sub- 
triangle may be further subdivided in like manner if 
not considered suitable for Gouraud shading, other- 
wise it is Gouraud shaded. Hence, surface normal 
vector interpolation, a Phong shading technique, is 

20 used to a limited extent, but when a sub-triangle is 
suitable for Gouraud shading intensity values are 
calculated at each vertex and the sub-triangle is 
tiled by linear interpolation of intensity values, i.e. a 
Gouraud shading technique. 

25 Triangle subdivision may be performed in 

world space or pixel space. When performed in 
world space the step of Gouraud shading must 
include a transformation function to pixel coordi- 
nates. When performed in pixel space the triangle 

30 vertices are converted to pixel coordinates prior to 
subdivision and the step of Gouraud shading is 
executed as previously described. In the following 
discussion it will be assumed that subdivision is 
performed in pixel space. 

as Because sub-triangles are formed by bisection, 

each edge of a sub-triangle is parallel to a cor- 
responding edge of its parent triangle. Accordingly, 
the values of DX/DY, s. Pine and Nine for each 
edge of a sub-triangle are the same as the values 

40 of DX/DY, s, Pine and Nine for a corresponding 
parallel edge of its parent triangle. Because values 
for these parameters are taken directly from the 
parent triangle, vertices of sub-triangles need not 
be converted to pixeJ coordinates to calculate 

45 Bresenham parameters and rounding errors asso- 
ciated with such conversions can be avoided. 

The sub-division algorithm is particularly useful 
as a recursive procedure whereby a parent triangle 
is the initial input parameter and its sub-triangies 

so are passed in subsequent recursive calls with re- 
curs ion being dependent upon a given sub-tri- 
angle's suitability for Gouraud shading. Thus, a 
parent triangle may be recursively subdivided into 
many smaller sub-triangles each corresponding to 

65 a smaller, more planar portion of the curved sur- 
face, and as such more suitable for Gouraud shad- 
ing than its parent Eventually, each sub-triangle is 
considered suitable for Gouraud shading and the 
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vertex data structures for its vertices must be com- 
pleted by calculating an AErr value (initial error 
value) and a seed grid point (Xo ( Y 0 ) for each sub- 
triangle edge to execute Gouraud shading. The 
seed grid point (Xo.Yo) indicates where tiling is to 
begin with respect to the associated edge and AErr 
indicates the distance in the X direction from 
(Xo.Yo) to the associated edge. 

To illustrate recursive subdivision of a parent 
triangle, consider triangle V0V1V2 as a parent tri- 
angle presented for tiling with each vertex of tri- 
angle V0V1 V2 defined in terms of positional coordi- 
nates and surface normal vectors. In the following 
discussion positional coordinate values are des- 
ignated with prefixes x, y and z while surface 
normal vector components are designated with pre- 
fixes a, b and c. Hence, the vertices of a parent 
triangle are associated with the following positional 
and surface normal vector parameters: 
V0 with (xO. yO. zO) and (aO, bO, cO) 
V1 with (x1, y1 t z1) and (a1, b1. d) 
V2 with (x2. y2, z2) and (a2, b2, c2) 

Prior to invoking the recursive procedure the 
vertices V0, V1. and V2 of the parent triangle are 
transformed to device (pixel) coordinates and the 
transformed vertices are sorted in ascending order 
according to the Y component. In the present ex- 
ample V0 has the smallest Y component and V2 
has the largest The Bresenham incremental pa- 
rameters for each edge of the parent triangle are 
computed, using the device coordinates of each 
vertex, and stored in the appropriate vertex data 
structure. 



For edge V1V0: 

OY10 = y1 - yO 

DX10 = x1 - xO 

slopelO = DX10/DY10 

slO = floor(slopelO) 

PlndO = DX10-s10*DY10- DY10 

NindO = DX10 -S10TDY10 

For edge V2V1: 

DY21 = y2 - y1 

DX21 = x2 - x1 

slope21 = DX21/OY21 

s21 = f!oor(slope21 ) 

Plnc21 = DX21 - s21 w DY21 - DY21 

Ninc21 = DX21 - s21*DY21 
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For edge V2V0: 

DY20 = y2 - yO 
DX20 = x2- xO 
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slope20 = DX20/DY20 

s20 = floor(slope20) 

Plnc20 = DX20 - s20"DY20 - DY20 

Ninc20 = DX20 - S20-DY20 

These Bresenham parameters are applicable to 
any sub-triangle of the parent, and the vertices of 
each sub-triangle are automatically sorted accord- 
ing to the Y component since the sub-triangles are 
similar to the parent triangle and the parent triangle 
is already sorted according to the Y component. 

The next step is to determine whether the 
parent triangle is suitable for Gouraud shading. As 
will be discussed in greater detail hereafter, several 
methods may be used to determine suitability for 
Gouraud shading. 

if the parent triangle is considered suitable for 
Gouraud shading, then the shading procedure dis- 
cussed is employed to shade the parent triangle. 
Otherwise, the vertex data structures of the parent 
triangle are passed to the following recursive pro- 
cedure entitled PSEUDO — PHONG. Local variables 
V01, V02, and V12 are vertex data structures cor- 
responding to the bisection points of edges V0V1 f 
V0V2 and V1V2. respectively. 

PSEUDO PHONG (VO. VI, V2) 

vertex V01, V02, V12 

begin pseudo__ phong 

t get bisection point V02 for edge V0V2 7 

V02.X : = (VO.x + V2.x)/2 

V02.y : = (VO.y + V2.y)/2 

V02.Z : = (VO.z + V2.z)/2 

A get surface normal vector for bisection point V02 
7 

V02.a:=(V0.a + V2.a) 

V02.b : = (VO.b + V2.b) 

V02.C : = (VO.c + V2,c) 

L : « NORM(V02.a, V02.b, V02.c) 

V02.a : = V02.a/L 

V02.b : = V02.b/L 

V02.C : = V02.C/L 

r get bisection point V01 for edge V0V1 7 
V01.X : = (V0.x + V1.x)/2 
V01.y := (VO.y + V1.y)/2 
V01.Z : = (V0.z + V1.z)/2 

t get surface normal vector for bisection point V01 
7 

V01 .a : =• (V0.a + V1 .a) 

V01.b := (VO.b + V1.b) 

VOI.c := (V0.c + V1.c) 

L:= NORM(V01.a. V01.b. V01.C) 

V01.a := V01.a/L 

V01.b := V01.b/L 

V01.C := V01.C/L 

t get bisection point V1 2 for edge V1 v*2 7 
V12.X : = (V1.x + V2.x)/2 
V12.y := (V1.y + V2.y)/2 
V12.z := (V1.z + V2.2)/2 

r get surface normal vector for bisection point V12 
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V12.a := (V1.a + V2.a) 

V12.b := (V1±> + V2.b) 

V12.c := (V1.c + V2.c) 

L := NORM(V12.a, V12.b t V12.C) 

V12.a := V12.a/L 

Vl2.b := V12.b/L 

V12.C := V12.C/L 

r test each sub-triangle to determine if suitable for 
Gouraud shading or if further subdivision is re- 
quired 7 

r for triangle V01V12V02 7 

if <RECURSE_CONDITION(V01 .V12.V02)) 

then PSEUDO_PHONG(V01,V12,V02) 

else GOURAUD SHADE (V01.V12.V02) 

rfor triangle V0V01V01 7 

if (RECURSE_CONDITION(V0,V01 ,V02)) 

then PSEUDO_PHONG(V0,V01,V02) 

else G O U R AU D SHAD E (V0.V01 ,V02) 

rfor triangle V01V1V12 7 

if <RECURSE_CONDITION(V01,V1.V12)) 

then PSEUDO — PHONG(V01,V1.V12) 

else GOURAUD_SHADE(V01.V1,V12) 

r for triangle V01 V12V2 7 

if <RECURSE_CONDITION(V02 f V1 2.V2)) 

then PSEUDO_ PHONG(V02,V12,V2) 

else GOURAUD SHADE (V01.V12.V2) 

end pseudo phong 

The depth of recursion, as determined by the 
boolean function RECURSE__CONDITION, may be 
established in an adaptive manner by examining 
the approximate curvature of the surface in the 
vicinity of the facet At each level of recursion, 
including the zeroth, a unit-length surface normal 
vector is given or calculated for each vertex of the 
triangle or sub-triangle. This is because unit-length 
normal vectors are needed for the lighting model 
computations. The three vertex surface normals of 
a triangle (or sub-triangle) are compared, by cal- 
culating their dot products and determining how 
close they are to unity. (The dot product of two unit 
vectors is equal to the cosine of the angle between 
the vectors.) If each of the three dot products 
exceeds a predetermined value, which depends on 
the degree of realism desired, It is assumed that 
the surface is planar in the vicinity of the triangle 
(or sub-triangle). Alternatively, the cross product of 
two vectors may be calculated and its magnitude 
compared to zero. When the magnitude of the 
cross product is sufficiently close to zero, the tri- 
angle or sub-triangle may be considered suitably 
planar and Gouraud shading of the triangle (or sub- 
triangle) swill provide a satisfactory rendering for a 
diffuse lighting model using point sources at infin- 
ity. 

In addition, recursion may be a function of the 
size of the sub-triangle in device coordinates. 
Since, at a given level of recursion, all sub-triangles 



are congruent, they are all of the same size and 
therefore use of a size limitation on recurs ion 
results in each sub-triangle being divided the same 
number of times. However, adjacent display plane 

s or parent triangles may be of different sizes, and 
therefore may not be divided the same number of 
times. When subdivision is performed in world 
space, this method introduces the possibility of 
gaps and overlaps between parent triangles as new 

70 vertices are introduced which should border adja- 
cent triangles, but do not. In one embodiment of 
the invention, the depth of recursion is fixed at four 
recursive calls for each display plane triangle, 
thereby dividing each display plane triangle into 

75 256 sub-triangles, and this has been found to pro- 
vide satisfactory results. 

A limitation on the depth of recursion is im- 
posed by the number of parameters that can be 
stored. Thus, on each subdivision additional param- 

20 eters are generated and must be stored for the 
next subdivision, and if the total number of param- 
eters that can be stored is limited, the number of 
subdivisions that can be performed is also limited. 
The function NORM returns the length or mag- 

25 nitude of a vector and may be implemented in 
several ways. A very accurate, but computationally 
expensive, version of NORM is as follows. 
NORM (a,b,c) = square_jroot(a*a + b # b + c*c) 
A faster, but slightly less accurate version is listed 

30 as 

NORM (a.b,c> = MAX (ABS(a),ABS(b).ABS(c)) 
where the function ABS returns the absolute value 
of its argument. Given the length L of a vector, a 
unit surface normal vector is obtained from the 

35 interpolated surface normal vector by dividing each 
component by L. 

If triangle subdivision is performed in world 

space then the procedure GOURAUD SHADE Is 

responsible for transforming triangle vertices into 

40 pixel coordinates. If subdivision is performed in 

pixel space, the procedure GOURAUD SHADE 

operates according to the previous description of 
Gouraud shading. 

Once a triangle or sub-triangle is determined to 

45 be suitable for Gouraud shading, the vertex data 
structures for its three vertices must be completed 
by calculating an AErr value and a seed grid point 
(Xo,Y 0 ) for each edge. Because vertices of a sub- 
triangle are not necessarily coincident with grid 

so points, calculation of these parameters is slightly 
different than the case of a parent triangle having 
its vertices at grid points. The Bresenham algo- 
rithm operates by incremental calculation of vari- 
ables to determine the position of a given grid 

55 point with respect to a line. In tiling a particular 
sub-triangle, it is necessary to determine values for 
AErr, Xo, and Y 0 as if the Bresenham algorithm 
were tiling the parent triangle sand had just 
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reached the particular sub-triangle. The Bresenham 
algorithm could be executed in simulation as if 
tiling the parent triangle in order to capture these 
values at the moment of reaching the region of the 
sub-triangle. However, the following calculation per- 
forms this task more efficiently. Consider edge 
V01V12 of triangle V01V12V02. From the vertex 
data structures for ths parent triangle, values for 
Pine, Nine, and s for edge V01V12 are known from 
the values for edge V0V2. Hence, 
slope = DX207DY20 
s = floor (slope) 

To calculate AErr and a seed grid point (X 0 .Y 0 ) 
for edge V01V12 to tile sub-triangle V01V12V02: 
Y 0 = f!oor(V01.y) 

Xoreal = V01.x - slope*(V01.y - YO) 
Xo - floor (Xoreal) 

AErr = floor (((Xo real + slope)-(Xo +s)-1)TDY20) 
If edge V01V02 is on the left, the starting X value 
for this edge is Xo + 1 , otherwise it is Xo. 

All necessary Bresenham parameters for edge 
V01 V12 are now available: s, Pine, and Nine from 
the parent triangle and AErr, Xo, and Yo as com- 
puted above. Similar calculations are performed to 
complete the Bresenham parameters for edges 
V12V02 and V02V01. Calculation of intensity values 
for vertices V01, V12, and V02 are derived by 
application of a selected lighting model to the re- 
spective vertex surface normal vectors and sub- 
triangle V01V12V02 is Gouraud shaded using XY 
address generator 2 and I generator 6 as described 
above. 

Thus a method of shading a display area has 
been shown which is less computationally expen- 
sive, yet nearly as good visually, compared to 
Phong shading. When used in conjunction with a 
Bresenham algorithm to shade sub-triangles 
formed by edge bisection, Bresenham incremental 
parameters calculated for the triangular area are 
applicable to its sub-triangles. 

It will be appreciated that the present invention 
is not restricted to the particular method that has 
been described, and that variations may be made 
therein without departing from the scope of the 
invention as defined in the appended claims and 
equivalents thereof. For example, it is not neces- 
sary that Gouraud shading be performed on a sub- 
triangle, since constant shading may be performed 
if a lower degree of realism is acceptable. In the 
case in which dot products are calculated in order 
to determine whether to make a further recursion, it 
is not necessary that all three dot products of the 
unit vertex normals for a triangle (or sub-triangle) 
be evaluated, since if two of the dot products are 
almost one, the third dot product will not be much 
smaller than one. 

The method of the present invention has ap- 
plication beyond triangle subdivision by edge 



bisection to obtain four congruent sub-triangles. 
Other methods may be used to generate sub- 
triangles. For example, parent triangles may be 
subdivided by determining the bisection point of 
5 one edge and connecting that point to the opposite 
vertex to divide a parent triangle into two sub- 
triangles, the positional coordinate values and sur- 
face normal vector component values being inter- 
polated from the endpoints of the bisected edge. 

70 

Claims 

1. A method of loading intensity values into a 
is memory device (5) having a plurality of addres- 
sable memory locations, each of which can be 
defined uniquely by an address having an X com- 
ponent and a Y component, which locations cor- 
respond respectively to points in a rectangular ar- 

20 ray, and wherein it is desired to load memory 
locations corresponding to points lying within an 
area defined by at least three vertices (X0.Y0), 
(X1.Y1) and (X2.Y2). said method comprising: 
defining points (X3.Y3) and (X4.Y4) on lines that 

25 join the point (X0.YO) with the points (X1.Y1) and 
(X2,Y2) respectively, such that a sub-triangle is 
defined by the points (X0.Y0), (X3.Y3) and (X4.Y4), 
calculating intensity values 10, 13 and 14 for the 
points (X0.Y0), (X3.Y3) and (X4,Y4) respectively, 

30 and 

loading the memory locations corresponding to the 
points within the sub-triangle with intensity values 
that are derived from at least one of 10, 13 and 14. 

2. A method according to claim 1, comprising 
as loading the memory locations corresponding to the 

points within the sub-triangle with intensity values 
that are derived from 10, 13 and 14 in combination. 

3. A method according to claim 2, comprising 
loading the memory locations corresponding to the 

40 points within the sub-triangle with intensity values 
that are derived from 10, I3 and 14 by linear inter- 
polation. 

4. A method according to claim 1, comprising 
defining the points (X3.Y3) and (X4,Y4) by bisect- 

45 ing the lines that join the point (X0.Y0) with the 
points (X1.Y1) and (X2.Y2). 

5. A method according to claim 1 , wherein said 
memory locations correspond to grid points of a 
display device. 

so 6. A method of loading intensity values into a 

memory device (5) having a plurality of addres- 
sable memory locations, each of which can be 
defined uniquely by an address having an X com- 
ponent and a Y component, which locations cor- 

55 respond respectively to points in a rectangular ar- 
ray, and wherein it is desired to load memory 
locations corresponding to points lying within a 
triangle defined by three vertices (X0.Y0), (X1,Y1) 
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and (X2.Y2), and to which normal vectors NO, N1 
and N2 are respectively assigned, said method 
comprising: 

defining points (X3.Y3) and (X4.Y4) on the lines 
that join the point (XO.YO) with the points (X1.Y1) 
and (X2.Y2) respectively, such that a sub-triangle is 
defined by the points (XO.YO) (X3.Y3) and (X4.Y4). 
deriving interpolated normal vectors N3 and N4 for 
the points (X3.Y3) and (X4.Y4) from the normal 
vectors NO, N1 and N2, and 

if the normal vectors NO. N3 and N4 bear a pre- 
determined relationship to each other, calculating 
intensity values !0, 13 and 14 for the points (XO.YO), 
(X3.Y3) and (X4.Y4) respectively and loading the 
memory locations corresponding to the points with- 
in the sub-triangle with intensity values that are 
derived from at least one of 10, 13 and (4. 

7. A method according to claim 6 t further com- 
prising, if the normal vectors NO, N3 and N4 do not 
bear said predetermined relationship to each other, 
defining points (X5.Y5) and (X6.Y6) on the lines 
that join the point (X0.Y0) with the points (X3.Y3) 
and (X4.Y4) respectively, such that a subtriangle is 
defined by the points (XO.YO). (X5.Y5) and (X6.Y6), 
deriving normal vectors N5 and N6 for the points 
(X5.Y5) and (X6.Y6) and. if the normal vectors NO, 
N5 and N6 bear a predetermined relationship to 
each other, calculating intensity values 10, 15 and 16 
for the points (X0.Y0). (X5.Y5) and (X6.Y6) respec- 
tively and loading the memory locations corre- 
sponding to the points within the sub-triangle de- 
fined by the points (XO.YO), (X5.Y5) and (X6.Y6) 
with intensity values that are derived from 10, 15 
and 16 by linear interpolation. 

8. A method according to claim 6, comprising 
determining whether the normal vectors NO, N3 
and N4 bear said predetermined relationship to 
each other by evaluating at least two dot products 
of the normal vectors NO, N3 and N4. 

9. A method according to claim 6, comprising 
determining whether the normal vectors NO. N3 
and N4 bear said predetermined relationship to 
each other by evaluating at least two cross pro- 
ducts of the normal vectors NO. N3 and N4. 

10. A method according to claim 6, wherein 
said memory locations correspond to grid points of 
a display device. 

1 1 . A method of shading a representation of a 
triangular surface having normal vectors associated 
with each vertex (VO. V1. V2) of the triangular 
surface, comprising: (a) defining three points 
(V01.V02.V1 2) each being on one of the three 
edges of the triangular surface, such that first, 
second, third and fourth sub-triangles 
(V0.V01 ,V02;V01 .VI ,V1 2;V02,V1 2,V2;V01 .V02.V1 2). 
are formed; each of the first, second and third sub- 
triangles having one of its three vertices in com- 
mon with a corresponding vertex of the triangular 



surface and two of said three points as its other 
two vertices, and the fourth sub-triangle having said 
three points as its vertices respectively, 

(b) deriving interpolated normal vectors for 
s said three points respectively from the normal vec- 
tors associated with the vertices of the triangular 
surface, and 

(c) if the interpolated normal vectors for the 
three vertices of a sub-triangle bear a predeter- 

10 mined relationship to each other, calculating inten- 
sity values for those three vertices respectively and 
utilizing the calculated intensity values to shade the 
sub-triangle. 

12. A method of shading a representation of a 
is curved surface modeled by a facet defined by 

three vertices (V0,V1 f V2), three straight lines join- 
ing those vertices, and normals (N0.N1.N2) to the 
curved surface associated with each vertex, com- 
prising: 

20 (a) if the normals bear a predetermined rela- 

tionship to each other, calculating intensity values 
for the three vertices respectively and utilizing the 
calculated intensity values to shade the facet, 
(b) otherwise: 

25 (i) calculating three points (V01.V02.V1 2) that bi- 
sect the three straight lines respectively, whereby 
four congruent sub-triangles 

(V0.V01 ,V02; V01 f V1 f V1 2;V02.V1 2. V2; V01 , V02.V1 2) 
are defined, 

30 (ii) -deriving interpolated normals (N01.N02.N12) for 
said three points respectively from the normals 
associated with each vertex defining the facet and 
(iii) if the interpolated normals for the three vertices 
of a sub-triangle bear a predetermined relationship 

as to each other, calculating intensity values for these 
three vertices respectively and utilizing the cal- 
culated intensity values to shade the sub-triangle. 

13. A method according to claim 12. wherein 
the method is recursive by application of the meth- 

40 od of claim 12 to shade a sub-triangle when the 
interpolated normals of the sub-triangle do not bear 
said predetermined relationship. 

14. A method of loading a memory device 
having a plurality of addressable memory locations, 

45 each of which can be defined uniquely by an 
address having an X component and a Y compo- 
nent, which locations correspond respectively to 
grid points in a rectangular array, and wherein it Is 
desired to load memory locations corresponding to 

so grid points lying within a display area defined by at 
least three vertices (XO.YO). (X1.Y1) and (X2.Y2), 
which are grid points, and corresponding surface 
normal vectors NO, N1 and N2, respectively, said 
method comprising: 

.55 defining a vertex (X3.Y3) by reference to at least 
two of vertices (XO.YO), (X1.Y1) and (X2.Y2), 
defining a surface normal vector N3 corresponding 
to vertex (X3.Y3) by interpolation of at least two of 



9 



17 EP 0 366 463 A2 18 



vectors NO. N1 and N2, 

defining a sub-area within said display area defined 
by at least three of vertices (XO.YO), (X1.Y1), 
(X2.Y2) and (X3.Y3). 

calculating intensity values for said vertices defin- 
ing said sub-area using the surface normal vectors 
corresponding to said vertices defining said sub- 
area, and 

loading the memory locations corresponding to grid 
points within said sub-area with intensity values 
that are derived from intensity values for said ver- 
tices defining said sub-area. 

15. A method according to claim 14 wherein 
said display area is a triangle defined by said 
vertices (XO.YO), (X1.Y1) and (X2.Y2), said vertex 
(X3.Y3) is positioned along a line connecting ver- 
tices (XO.YO) and (X1.Y1) and derived by interpola- 
tion of vertices (XO.YO) and (X1.Y1), the surface 
normal vector corresponding to vertex (X3.Y3) is 
derived by interpolation of the surface normal vec- 
tors corresponding to vertices (XO.YO) and (X1.Y1), 
and the method further comprises the steps: 
defining a vertex (X4.Y4) positioned along a line 
connecting vertices (XO.YO) and (X2.Y2) by inter- 
polation of vertices (XO.YO) and (X2.Y2), 
deriving a surface normal vector corresponding to 
vertex (X4.Y4) by interpolation of the surface nor- 
mal vectors corresponding to vertices (XO.YO) and 
(X2.Y2), 

calculating an intensity value for vertex (X4.Y4) 
using the surface normal vector corresponding to 
vertex (X4.Y4) , and whereby said sub-area is a 
triangle defined by vertices (XO.YO), (X3.Y3) and 
(X4.Y4) and intensity values corresponding to grid 
points within said sub-triangle are derived by inter- 
polation of intensity values for vertices (XO.YO), 
(X3.Y3) and (X4.Y4). 

16. A method according to claim 14 wherein 
said display area is a triangle defined by said 
vertices (XO.YO), (X1.Y1) and (X2.Y2). said vertex 
(X3.Y3) is positioned along a line connecting ver- 
tices (XO.YO) and (X1.Y1) and derived by interpola- 
tion of vertices (XO.YO) and (X1.Y1), the surface 
normal vector corresponding to vertex (X3.Y3) is 
derived by interpolation of the surface normal vec- 
tors corresponding to vertices (X0,Y0) and (X1.Y1), 
said sub-area is a triangle defined by vertices 
(XO.YO). (X2.Y2) and (X3.Y3), and intensity values 
corresponding to grid points within said sub-tri- 
angle are derived by interpolation of the intensity 
values for vertices (XO.YO), (X2.Y2) and (X3.Y3). 

1 7. A computer graphics system comprising: 

(a) a display device (7) having a display 
surface. 

(b) display address means (8) for generating 
a display address signal having an X component 
and a Y component such as to define a set of grid 
points on the display surface in a rectangular array, 



and 

(c) picture processor means (1.6) for receiv- 
ing a display list containing the locations of at least 
three vertices (XO.YO). (X1.Y1) and (X2.Y2). the 
s picture processor means automatically generating 
intensity values corresponding to points lying within 
an area defined by the three vertices by performing 
the following operations, 

(i) defining points (X3.Y3) and (X4.Y4) on lines that 
w join the point (XO.YO) with sthe points (X1.Y1) and 

(X2.Y2) respectively, such that a sub-triangle is 
defined by the points (XO.YO), (X3.Y3) and (X4.Y4), 

(ii) calculating intensity values 10, 13 and 14 for the 
points (XO.YO), (X3.Y3) and (X4.Y4) respectively, 

is and 

(in) calculating intensity values for points within the 
sub-triangle from at least one of 10, 13 and 14. 

18. A computer graphics system according to 
claim 17, comprising an electronic memory device 

20 (5) having a set of addressable memory locations 
equivalent to the set of grid points and correspond- 
ing respectively to the grid points, and means 
(2,3,4) for loading the intensity values calculated for 
points within the sub-triangle into the correspond- 

25 ing memory locations of the memory device. 

19. A computer graphics system according to 
claim 18, wherein the means for loading intensity 
values into the memory device comprise an XY 
address generator (2) connected to the picture 

30 processor means for automatically generating the 
addresses of the memory locations corresponding 
to points within the sub-triangle and applying the 
addresses to the memory device. 

20. A computer graphics system according to 
35 claim 17, wherein the picture processor means 

comprise a picture processor (1) for receiving the 
data structures and an I generator (b) connected to 
the picture processor for generating the intensity 
values. 

40 21. A computer graphics system according to 

claim 17, wherein the picture processor means 
calculate intensity values for points within the sub- 
triangle from 10, 13 and 14 in combination. 

22. A computer graphics system according to 
45 claim 21, wherein the picture processor means 

calculate intensity values for points within the sub- 
triangle from I0 ( 13 and 14 by linear interpolation. 

23. A computer graphics system according to 
claim 17, wherein the picture processor defines the 

so points (X3.Y3) and (X4.Y4) by bisecting the lines 
that join the point (X0.Y0) with the points (X1.Y1) 
and (X2.Y2). 

24. A computer graphics system comprising: 
(a) a display device (7) having a display 

55 surface, (b) display address means (8) for generat- 
ing a display address signal having an X compo- 
nent and a Y component such as to define a set of 
grid points on the display surface in a rectangular 
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array, and 

(c) picture processor means (1 .6) for receiv- 
ing a display list containing the locations of at least 
three vertices (XO.YO). (X1.Y1) and (X2.Y2) and 
normal vectors NO, N1 and N2 for the three ver- 
tices (XO.YO), (X1.Y1) and (X2.Y2) respectively, the 
picture processor means automatically generating 
intensity values corresponding to points lying within 
an area defined by the three vertices by performing 
the following operations, 

(i) defining points (X3.Y3) and (X4,Y4) on lines that 
join the point (XO.YO) with the points (X1.Y1) and 
(X2.Y2) respectively, such that a sub-triangle is 
defined by the points (XO,YO), (X3.Y3) and (X4.Y4), 

(ii) deriving normal vectors N3 and N4 for the 
points (X3,Y3) and (X4.Y4) from the normal vectors 
NO, N1 and N2, 

(iii) determining whether the normal vectors NO, N3 
and N4 bear a predetermined relationship to each 
other and, if so, 

(iv) calculating intensity values 10, 13 and 14 for the 
points (X0,Y0) . (X3.Y3) and (X4,Y4) respectively, 
and 

(v) calculating intensity values for points within the 
sub-triangle from at least one of 10. 13 and 14. 

25. A computer graphics system according to 
claim 24, comprising an electronic memory devices 
(5) having a set of addressable memory locations 
equivalent to the set of grid points and correspond- 
ing respectively to the grid points, and means 
(2,3,4) for loading the intensity values calculated for 
points within the sub-triangle into the correspond- 
ing memory locations of the memory device. 

26. A computer graphics system according to 
claim 25, wherein the means for loading intensity 
values into the memory device comprise an XY 
address generator (2) connected to the picture 
processor means for automatically generating the 
addresses of the memory locations corresponding 
to points within the sub-triangle and applying the 
addresses to the memory device. 

27. A computer graphics system according to 
claim 24, wherein the picture processor means 
comprise a picture processor (1) for receiving the 
data structures and an I generator (6) connected to 
the picture processor for generating the intensity 
values. 

28. A computer graphics system according to 
claim 24, wherein the picture processor means 
determine whether the normal vectors NO, N3 and 
N4 bear said predetermined relationship to each 
other by evaluating at least two dot products of the 
normal vectors NO, N3 and N4. 

29. A computer graphics system according to 
claim 24, wherein the picture processor means 
determine whether the normal vectors NO, N3 and 
N4 bear said predetermined relationship to each 
other by evaluating at least two cross products of 
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the normal vectors NO, N3 and N4. 

30. A computer graphics system comprising: 
(a) a display device (7) having a display 
surface, 

s (b) display address means (8) for generating 

a display address signal having an X component 
and a Y component such as to define a set of grid 
points on the display surface in a rectangular array, 
and 

io (c) picture processor means (1 ,6) for receiv- 

ing a display list containing the locations of three 
vertices defining a triangular surface and three nor- 
mal vectors associated with the vertices respec- 
tively, the picture processor means shading a re- 

is presentation of the triangular surface by performing 
the following operations, 

(i) defining three points each being on one of the 
three edges of the triangular surface, such that 
three sub-triangles are formed, each sub-triangle 

20 having one of its three vertices in common with a 
corresponding vertex of the triangular surface and 
two of said three points as its other two vertices, 
and a fourth sub-triangle is formed with said three 
points as its vertices respectively, 

25 (ii) deriving interpolated normal vectors for said 
three points respectively from the normal vectors 
associated with the vertices of the triangular sur- 
face, and 

(iii) if the interpolated normal vectors for the three 
30 vertices of a sub-triangle bear a predetermined 
relationship to each other, calculating intensity val- 
ues for those three vertices respectively and utiliz- 
ing the calculated intensity values to shade the 
sub-triangle. 

3s 31. A computer graphics system according to 

claim 30, comprising an electronic memory device 
(5) having a set of addressable memory locations 
equivalent to the set of grid points and correspond- 
ing respectively to the grid points, and means 

40 (2,3,4) for loading the intensity values calculated for 
points within the sub-triangle into the correspond- 
ing memory locations of the memory device. 
32. A computer graphics system comprising: 

(a) a display device (7) having a display 
45 surface, 

(b) display address means (8) for generating 
a display address signal having an X component 
and a Y component such as to define a set of grid 
points on the display surface in a rectangular array, 

so and 

(c) picture processor means (1,6) for receiv- 
ing a display list containing the locations of at least 
three vertices (XO.YO) . (X1.Y1) and (X2.Y2) defin- 
ing a display area, and surface normal vectors NO, 

55 N1 and N2 associated respectively with the ver- 
tices (XO.YO). (X1.Y1) and (X2.Y2), the picture pro- 
cessor means automatically generating intensity 
values corresponding to points lying within an area 
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defined by the three vertices by performing the 
following operations, 

(i) defining a vertex (X3.Y3) by reference to at least 
two of vertices (XO.YO), (X1.Y1) and (X2.Y2), 

(ii) defining a surface normal vector N3 correspond- s 
ing to vertex (X3.Y3) by interpolation of at least two 

of vectors NO. N1 and N2, 

(iii) defining a sub-area within said display area 
defined by at least three of vertices (XO.YO), 

(X1 t Y1 ) f (X2.Y2) and (X3 ( Y3) t 1Q 

(iv) calculating intensity values for said vertices 
defining said sub-area using the surface normal 
vectors corresponding to said vertices defining said 
sub-area, and 

(v) calculating intensity values for grid points within rs 
said sub-area from the intensity values for said 
vertices. 

33. A computer graphics system according to 
claim 32, comprising an electronic memory device 
(5) having a set of addressable memory locations 20 
equivalent to the set of grid points and correspond- 
ing respectively to the grid points, and means 
(2,3,4) for loading the intensity values calculated for 
points within the sub-triangle into the correspond- 
ing memory locations of the memory device. 25 
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